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Abstract 

■ We propose a new point of view for interpreting Newton's and Ein- 
stein's theories of gravity. By taking inspiration from Continuum Me- 

^f*j ] chanics and its treatment of anisotropics, we formulate new theories of 

£N| ■ gravity- These models are simple and natural generalisations with many 

interesting properties. Above all, their precise form can, in principle, be 
determined experimentally. 

^! 1 Introduction 

Since Newton's formulation of the law of universal gravitation in the Principia 
1687, the theory has been unchanged. Very few attempts have been made to 
modify Newtonian gravity, modified Newtonian dynamics being probably the 
. one exception pQ. It took almost 330 years until a generalisation of Newton's 

■^J- | theory was successfully constructed. Einstein's theory of gravity was radically 

i different to any previous physical theory. It abandoned the absolute notion of 

time and replaced the concept of force by the curvature of a four dimensional 
space. However, it took only a few years when the first modifications and ex- 
tensions of Einstein's theory appeared in the literature, see [2] for an excellent 
historical review. Ever since, modified theories of gravity have enjoyed a promi- 
nent role in theoretical physics. It is the aim of the current work to propose 
modifications of Newton's and Einstein's theories of gravity by applying the 
same idea to both of them. 

■ Let us start with the actions of Newtonian gravity and its relativistic ana- 
logue, the Einstein-Hilbert action 

(1) 
(2) 

Variation of the Newtonian action with respect to the gravitation potential (p 
yields the well-known Poisson equation Atp = AirGp where p is a given matter 
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distribution. Likewise, variations of the Einstein-Hilbert action with respect 
to the metric (which contains the gravitational potentials) gives the famous 
Einstein field equations G a p — 87rG/c 4 T Q( g. 

When comparing the Einstein-Hilbert action (0) with other models in physics, 
it appears to be somewhat unnatural as one generally considers potential ener- 
gies quadratic in the field strength, Hooke's law probably being the best known. 
However, when looking at the Einstein-Hilbert action, there are various very 
good reason for its form. Some of these can be motivated mathematically while 
others are simply observational. The field equations derived from it are close 
to Newton's theory of gravity and where the solutions deviate they are doing 
so in precisely the way to be in agreement with observations. However, there 
are observational facts which strongly indicate that our understanding of the 
gravitational force is far from complete. Dark matter and dark energy are two 
unknown forms of matter which are required to make the Universe work. 

Looking at (jTJ) and © we note that both actions contain a contraction 
with respect to the metric. In the Newtonian model we simply have the flat 
metric for Cartesian 3-space. We are now taking some inspiration from Con- 
tinuum Mechanics and in particular constitutive equations which specify the 
various different models studied in this field. In Maxwell's electromagnetism 
the constitutive equations for instance define the form of the dielectric ten- 
sor D = eE where e is a rank 2 tensor in general and we should really write 
Di = e^Ej. When working with the Faraday tensor and its corresponding 
excitation, one writes Hij — 1/4 £ij m nX mnkl Fki where for Maxwell's theory in 
vacuum ~^ mnkl = y/—g(g mk g nl — g nk g ml ), see [3J- In many simple applications 
£j? is taken to be proportional to the Kronecker delta which corresponds to 
a simple isotropic medium with dielectric constant e, namely — edf . In 
other words, a rank 2 isotropic tensor is proportional to the metric tensor. We 
are using this observation to argue that the Newton and Einstein-Hilbert ac- 
tions are based on the assumption of an isotropic medium. Since Nature has 
a strong tendency to be anisotropic, this appears to be rather unnatural. By 
doing so, we are changing the interpretation of the terms 5^ and g^ v in (fTJ) 
and ([2]), respectively. We now view them much like the material tensors speci- 
fied by constitutive equations which define the properties of the material. It is 
an interesting historical fact that F. Klein noted in a letter to Hilbert, as early 
as 1917, that the Ricci scalar in the Einstein-Hilbert action can be written as 
X mnkl Rmnkii see [1]. For a recent paper inspired by similar thoughts, see [5]. 

We might want to, for sake of concreteness, speak of the properties of the 
vacuum here. However, we would like to be very careful and point out that this 
change of viewpoint has a variety of philosophical implications when interpreting 
the new theories. We would like to keep these issues aside for now and proceed 
with the formulation of the theories. The assumption of an isotropic vacuum is 
certainly well supported by a host of experimental observations. However, there 
is clearly room for some improvement, in particular in general relativity. Many 
of the well-known modifications are is some sense rather severe. Field equations 
become higher than second order, new fields are introduced, locality is broken, 
local Lorentz invariance is broken etc. We will show that our new theories are 
completely harmless and retain all desired properties. 
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2 New theories of gravity 



One of the most conservative modifications of Newton's and Einstein's theories 
we can think of is therefore the following 



*~*Newton 



-p^+^-r, c ' 3 9 i <P d i <p\d 3 x, (3) 



8ttG~ " ,r 3 . 



9d"x, (4) 



where C lJ and are symmetric rank 2 tensors which contain information 
about the underlying structure of the theory. In Continuum Mechanics such 
objects are often referred to as material tensors or elastic coefficients [B|. We will 
stick with this well established notation and should clearly distinguish it from 
the energy-momentum tensor T a p. We should emphasise that in this approach 
to gravity we are able to assume any symmetry for the metric, for instance 
spherical symmetry or homogeneity and isotropy. However, the symmetry of 
the metric is independent of the symmetry of the elastic coefficients. 

Let us firstly consider variations of the new Newtonian theory §5§ with re- 
spect to the potential <f>. We find 

8 i (C ii d J <p)=4!rGp, (5) 

which reduces to the Poisson equation if we choose C ,J = 5^ . As C 1 -? can in 
principle be an arbitrary tensor, solutions to this equation may be quite different. 
As a simple example, let us consider the case when C 4 -? = diag(ci, C2, C3) with 
the Ci being constants. In continuum mechanics or solid state physics one would 
speak of a crystal with three principal propagation directions. In this case, the 
new field equation §5§ can be reduced to the Poisson equation by rescaling the 
coordinates. Looking for a radially symmetric solution gives the interesting 
result 

ip oc = . (6) 

V c\x 2 + c 2 y 2 + c 3 z 2 

In this case the gravitational field of a massive body is ellipsoidal instead of 
being spherical. However, the strength of this effect depends on the values of the 
constants. If for instance, the numerical values of the three constants ci, 02,03 
are very close to each other, then the gravitational field will look spherical unless 
very large distances are taking into account. 

It should also be noted that the components of C IJ do not have to be con- 
stants. They can be functions of the coordinates. Thus, the form of the gravi- 
tational field may be different in different regions of space, however, the grav- 
itational law ([5]) would still be universal. If we consider C 13 = x(t, x, y, z)5 %: > , 
the new field equations are 

Acp= p . (7 

X X 

If the function \ is slowly varying in space \dix\ *C I, we would have a the- 
ory which would be in very good agreement with Newton's theory on smaller 
scales, like the solar system. Note also that assuming a time varying x would 
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correspond to a time dependent gravitational constant in agreement with Dirac 
large number hypothesis [7]. 

Let us now consider the new action Q . We assume that the tensor C^ v may 
depend on the metric, this is important so that the limit of general relativity 
can be recovered. Note that the variation of the Ricci tensor does not give 
a surface term anymore because C^ v is not the metric. One has to integrate 
by parts twice (the Ricci tensor contains second derivatives of the metric) and 
move the derivatives action on the metric variations to act on the matter tensor. 
Denoting 

iff"' 
og a /3 

the whole calculation yields the field equation 

+ \g a ^^ v c^ v - V M V (Q C' 3 ^ = ^f-T afs . (9) 

where we have also included the energy- momentum tensor of matter. 

Recall that in General Relativity the field equations imply the conservation 
equations by virtue of the twice contracted Bianchi identities. This does no 
longer hold due to the presence of the material tensor in the field equations. 
Nonetheless, we can take the covariant derivative of the field equation ^ with 
respect to V a and find the following conservation equation 



- C aa V a R a - R <7 V a C a ' 7 = ^V a T a0 . (f 0) 

c 4 

This implies that the conservation equations for T Q/9 is no longer a consequence 
of the gravitational field equations and some additional information is required 
to close the system of equations. However, in the case when a vacuum is con- 
sidered, T a P — 0, the conservation equation (fTUf imposes constraints on the 
components of C^ v and thus its components cannot be specified arbitrarily. 

It should also be emphasised that one has to be careful when choosing the 
tensor and the derived quantity £^ q p , Namely, there is a conceptual dif- 
ference between prescribing the tensor C£ from prescribing C^ u . In the former 
case, we find that C^ v = g^C^ and thus C^ v has an explicit dependence on 
the metric which in turn will yield a non-trivial form for Yi^ va ^ . On the other, 
if we specify a priori, then this tensor does not depend on the metric and 
hence the resulting variational derivative would be zero. One has to carefully 
distinguish both cases when analysing the field equations. 

It is very difficult to analyse the field equations © given general S^"/ 3 
and C^ v . It is also not clear what would constitute a good choice for those 
quantities. As a first attempt to understand the theory, we assume & v to be 
conformally related to the inverse metric 

= <j)(x a )g^, (II) 
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where <j>(x a ) is a scalar function depending on the coordinates. This choice 
simplifies the field equation considerably and they are now given by 

4>G a0 +g afi U4>~ V a V p <t> = ^±T aP . (12) 

c 4 

Upon division by the field <p and solving for the Einstein tensor, the resulting 
equation 

Qap = 8^ TQ/3 + lya V Pf _ g*PU<j>) . (13) 

shows similarities with non-minimally coupled scalar field theories see Eq. (11)] 
The main difference between the current theory and most other approaches is 
that <j) is not a dynamical degree of freedom because C^ v is a prescribed tensor 
and has not corresponding equations of motion. The conformal model appears 
to agree with [5] in the limit when the Brans-Dicke parameter uo — > 0. How- 
ever, as variations with respect to <f> are not considered, there is no propagation 
equation for the scalar field. By substituting (fTTj) back into the action ((4]), we 
note that this simply corresponds to assuming the gravitational constant to be 
varying in time and space, see for instance [7J. Varying constants models are 
generally based on non-minimally coupled scalar field with kinetic term similar 
to Brans-Dicke theory, see 



3 A Schwarzschild like solution 

Let us start by considering a static and spherically symmetric vacuum spacetime 
described by the metric 

ds 2 = -e u ^dt 2 + e^dr 2 + rW , (14) 

where dfl 2 is the usual line element of the two-sphere. We also assume <f> — e^ r ). 
When the analogue situation is analysed in General Relativity where £ = 0, one 
finds two independent equations which determine the two unknown functions 
v(r) and /i(r). In this model, there is the additional degree of freedom £ and 
fortunately, there are now three independent equations. These are given by 



r i 2 r r 



e M 1 1 v 1 £' 

r z r L r r 

+f + e" + (e') 2 + - = o. a?) 

Zr r 

Eqn. (1161) can be solved for the function /i which can then be substituted into 
the other two equations. Combining those linearly gives the condition v 1 oc £' 
which then allows us to reduce this problem to a single differential equation. 
While separation of variable and subsequent integration is possible, the resulting 
equation cannot be solved analytically for the unknown functions due to its high 
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nonlincarity. However, by assuming £' = —C 2 v' with C < 1 we can find an 
approximate solution to the field equations which is given by 



e v = 1 - - + (9(C 4 ) , (18) 
r 

C 9C 3 

e-" = l-- + — + 0(C 4 ), (19) 
r r 

e £ = 1 + — + 0(C 4 ) . (20) 
r 

One easily verifies that this solution satisfies the field equations (|15 [l — ([TT |) up 
to 0(C A ). By choosing C = 2GM we arrive at a Schwarzschild like solution 
with only a small difference. Clearly, this difference of the order of C 3 and 
therefore it would be very difficult to distinguish between this metric and the 
Schwarzschild metric using solar system tests. This is a promising result which 
indicates that this theory can pass solar system tests without great difficulty. 



4 Cosmology 

Next, we want to study the cosmological implications of field equations 
Similar to the above we assume C^ u = <j){t) g^" and consider a flat FLRW 
universe 

ds 2 = -dt 2 + a(t)(dx 2 + dy 2 + dz 2 ) , (21) 

where a(t) is the scale factor. This yields the following cosmological field equa- 
tions 

ZH 2 cj> + ZH^=^- P , (22) 
c 4 

■■a n SttG 

-(t>-24>--2H4>-H 2 cl) = — r p, (23) 
a c 4 

where an overdot denotes differentiation with respect to time derivative and H 
is the Hubble parameter H = a/ a. Moreover, the conservation equation (|10[) 
implies 

^f{p + 3H(p + p)}=34>H + 6H 2 <j ): (24) 
c 4 

which can also be derived from the two field equations. 

We start looking for vacuum solutions, p = p = 0. Eqn. (f2"2"j) yields H = 
—(j>/4> which implies 

</>cxi. (25) 

Next consider (l23l) . inserting relation (|25|) gives the simple equation aa+(a) 2 = 
which is solved by a(f) = a^t — to which coincides with the evolution of a 
radiation dominated universe in standard cosmology. This means that in this 
simple model the vacuum energy naturally leads the universe into a radiation 
dominated era. 
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Under the assumption that C 7 ' 1 ' = 4>{t) g^ u and p = p = we cannot find 
an accelerating solution to the field equations. We could, of course, add reg- 
ular matter to the field equations and seek other solutions. However, we are 
particularly interested in the vacuum equations without ordinary matter. The 
reason for this is simply that we want to show that an additional non-dynamical 
structure in the theory suffices to get a dynamic universe. One can easily find 
such solutions by introducing an extra degree of freedom in the material tensor. 
Let us choose the elastic coefficients to be 

C£ = diag(-e,o-,a-,(7), (26) 

which means that C^ v — g^ a C^ and thus we have an explicit dependence on the 
metric, implying a non-zero given by Y^ va $ = g m C^ . One can interpret 

the quantities g and a as the energy density and pressure of the vacuum, thereby 
specifying its internal structure. 

The cosmological field equations of this model are 

2f4 + (<?-<7)-=0, (27) 
cr a 

a 2 a a . . 

2g— + 4g- + 4g- + {g - a) = , (28 
a z a a 

where the conservation equations can be derived from these two independent 
equations. In order to close the system, we choose a linear equation of state 
cr = wg. Now, Eq. (f2"T]) can be integrated and yields 

g = Ca 2/{w - 1 \ (29) 

where C is some constant of integration. We can substitute this into ((28)) and 
find 

a = a (t-t )r^ . (30) 

By virtue of (|29p we arrive at 

{t-t y^. (31) 

These are expanding and accelerating solutions provided < w < 1/3. The 
case w — 1/3 needs to be treated separately as the exponents tend to become 
very large. One would expect those to correspond to exponential functions, and 
indeed in this case 

a = a e xt , g = g Q e- 3Xt , (32) 

where go = A/oq. Thus, we were able to find solutions of the field equations 
modelling a universe which can undergo periods of accelerated expansion, with- 
out the need to introduce any forms of matter. All we have done is to add an 
additional non-dynamical structure to the theory on a very fundamental level. 

5 Conclusions 

It is tempting to argue that we introduced a form of matter through the back 
door by choosing our material tensor, the elastic coefficients. However, it is far 
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from clear whether this is indeed the case. Note that C^ v is not a dynamical 
variable and thus it cannot be interpreted simply as matter. In Continuum 
Mechanics and when working with crystal symmetries, the tensor C^ v is said 
to encode the symmetry properties of the material, in our case the vacuum. 
We simply say that the vacuum as we know it may have an internal structure 
which is specified by C^ v . We are breaking away from the assumption that the 
vacuum is isotropic and structureless. 

The form of the elastic coefficients can in principle be determined observa- 
tionally. In the context of cosmology, one could start with the specific given 
by ((26|) and assume it to be close to the metric g^ v . When one considers models 
where deviations from General Relativity will vary with cosmological time, it 
would be most interesting to see how observational data would determine the 
form of C^" which provides the best fit to the data. Using our approach to 
gravity, we will be able to use observations directly to specify the model instead 
of guessing new theories and deriving their implications on a case by case basis. 
Ultimately, experiments and observations will determine the correct theory. 
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